<!DOCTYPE html>
<html>

<head>
  <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
  <meta name="viewport"
    content="width=device-width,initial-scale=1,user-scalable=0,minimum-scale=1.0,maximum-scale=1.0" />
  <meta name="author" content="火星科技 http://mars3d.cn ">
  <meta name="apple-touch-fullscreen" content="yes">
  <meta name="apple-mobile-web-app-capable" content="yes">
  <meta name="apple-mobile-web-app-status-bar-style" content="black">
  <meta name="format-detection" content="telephone=no">
  <meta name="x5-fullscreen" content="true">
  <meta http-equiv="X-UA-Compatible" content="IE=Edge,chrome=1" />
  <!-- 标题及搜索关键字 -->
  <meta name="keywords" content="火星科技,cesium,3D,GIS,marsgis,三维,地球,地图,开发,框架,系统,示例,资料,模型,离线,外包,合肥,安徽,中国" />
  <meta name="description"
    content="火星科技 合肥火星 合肥火星科技 合肥火星科技有限公司 leaflet leaflet框架 leaflet开发 cesium cesium开发 cesium框架 三维 地球 模型  gis marsgis 地图离线 地图开发 地图框架 地图外包 框架 开发 外包  地图离线 二维地图 三维地图 全景漫游 地理信息系统 云GIS 三维GIS GIS平台 WebGIS" />

  <link rel="shortcut icon" type="image/x-icon" href="http://mars3d.cn/favicon.ico">
  <title>三维地图 | Mars3D | 三维地图 | 火星科技 | 合肥火星科技有限公司</title>


  <!--第三方lib-->
  <script type="text/javascript" src="../lib/include-lib.js" libpath="../lib/"
    include="jquery,font-awesome,bootstrap,layer,haoutil,turf,mars3d,kriging"></script>
  <link href="css/style.css" rel="stylesheet" />
</head>

<body class="dark">
  <!--加载前进行操作提示，优化用户体验-->
  <div id="mask" class="signmask" onclick="removeMask()"></div>


  <div id="mars3dContainer" class="mars3d-container"></div>


  <script src="js/kriging_config.js"></script>

  <script src="./js/common.js"></script>
  <script type="text/javascript">
    'use script' //开发环境建议开启严格模式

    var map

    function initMap(options) {
      //合并属性参数，可覆盖config.json中的对应配置
      var mapOptions = mars3d.Util.merge(options, {
        scene: {
          center: { lat: 23.057434, lng: 110.106285, alt: 8580407, heading: 356, pitch: -82 },
        },
      })

      //创建三维地球场景
      map = new mars3d.Map('mars3dContainer', mapOptions)

      $.getJSON('http://data.marsgis.cn/file/geojson/temperature.json', function (data) {
        //获取插值绘制图片
        var image = loadkriging(data.features, kriging_bounds, kriging_colors) //kriging_bounds, kriging_colors在js/temperature_data.js中定义

        var tileLayer = new mars3d.layer.ImageLayer({
          url: image,
          rectangle: {
            xmin: 73.4766,
            xmax: 135.088,
            ymin: 18.1055,
            ymax: 53.5693,
          },
          alpha: 0.4,
        })
        map.addLayer(tileLayer)
      })
    }

    //根据 克里金法 插值绘制图片
    function loadkriging(tempture, bounds, colors) {
      var canvas = document.createElement('canvas')
      canvas.width = 2000
      canvas.height = 2000

      var t = []
      var x = []
      var y = []
      for (var i = 0, len = tempture.length; i < len; i++) {
        t.push(tempture[i].properties.Temperatur) // 权重值
        x.push(tempture[i].geometry.coordinates[0]) // x
        y.push(tempture[i].geometry.coordinates[1]) // y
      }
      //1.用克里金训练一个variogram对象
      var variogram = kriging.train(t, x, y, 'exponential', 0, 100)

      //2.使用刚才的variogram对象使polygons描述的地理位置内的格网元素具备不一样的预测值；
      //bounds:普通的geojson格式的面的格式的coordinates。
      var grid = kriging.grid(bounds, variogram, 0.05)
      //3.将得到的格网预测值渲染到canvas画布上
      kriging.plot(canvas, grid, [73.4766, 135.088], [18.1055, 53.5693], colors)

      var image = canvas.toDataURL('image/png')
      return image
    }
  </script>
</body>

</html>
